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NOMEN  CU  TURK 


^  Dimensionless  distance  from  the  midchord  aft  to  the  axis  of 

rotation 

b  Semi-chord  length  of  foil 

C  Damping  coefficient 

Critical  damping  coefficient 
C(k)  Theodorsen's  function  of  reduced  frequency 

f  Carriage  excitation  frequency 

g  Artificial  damping  coefficient 

h  Translational  displacement  of  the  axis  of  rotation  (see  Figure 

4).  When  used  as  a  subscript,  h  refers  to  the  translational 
degree  of  freedom 

h^  Amplitude  of  the  oscillatory  translational  displacement 

I  Mass  moment  of  inertia  of  the  rotational  assembly  with  respect 

to  its  axis 


K  Linearized  spring  constant  of  the  rotational  degree  of  freedom 

CL 

K,  Linearized  spring  constant  of  the  translational  degree  of  freedom 
n 

k  Reduced  frequency  of  vibration,  k  =  tJUb/U 

L  Unsteady  lift 

L,  Dimensionless,  complex,  unsteady  lift  coefficients 

h  j  d 

i  Span  length  of  the  hydrofoil 

M  Unsteady  hydrodynamic  moment  about  the  axis  of  rotation 


\,a 


Dimensionless,  complex,  unsteady  moment  coefficients 

Mass  of  the  system  which  vibrates  in  translation  (including  the 
rotational  mass) 


r 


a 


Dimensionless  radius  of  gyration  about  the  elastic  axis 
r^  =  (l/rab^)^ 
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S  First  moment  of  mass  aft  of  the  rotational  axis 

a 

U  Velocity  of  the  foil 

X  Dimensionless  distance  aft  from  the  axis  tD  the  center  of  gravity 

of  the  rotational  assembly:  x  =  S  /mb 

a  a 

Z  Complex  parameter  used  in  flutter  determinant,  Z  =  (o)  /u))^(l  +  ig) 

CL 

a  The  rotational  displacement  of  the  foil  about  the  axis  of  rota¬ 

tion.  When  used  as  a  subscript,  it  refers  to  the  rotational  de¬ 
gree  of  freedom  (see  Figure  4) 

Amplitude  of  the  oscillatory  rotational  displacement 
p,  Mass-density  ratio;  p  =  m/TTpb^ i, 

p  Density  of  the  liquid 

Q  Ratio  of  the  translational  frequency  to  the  rotational  frequency, 

n  =  01,  /(ju 
h  a 

0)  Circular  frequency  of  the  ’coupled  vibration  in  water 

'0)^  Circular  frequency  of  vibration  in  the  rotational  mode  in  air 

OJ^  Circular  frequency  of  vibration  in  the  translational  mode  in  air 
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ABSTRACT 


Initial  evidence  of  flutter  of  a  fully  submerged 
hydrofoil  of  small  aspect  ratio  under  controlled  experi¬ 
mental  conditions  is  presented.  The  influence  of  several 
primary  parameters  on  flutter  was  investigated.  Theodorsen's 
two-dimensional,  unsteady  flutter  theory  yielded  flutter 
speed  predictions  in  good  agreement  with  experimental 
data.  Based  on  the  findings,  some  conditions  under  which 
flutter  of  displacement  or  hydrofoil  craft  could  be 
anticipated  are  discussed. 


INTRODUCTION 

Since  flutter  is  generally  a  destructive  phenomenon,  there  is  some 
concern  about  its  possible  occurrence  on  control-surfaces  in  water.  This 
concern  is  warranted  by  the  increasing  use  of  hydrofoil  configurations  at 
high  speed.  Several  Investigators  have  expressed  their  views  on  this 
matter  along  with  results  of  analytical  studies  in  the  new  field  of  hydro¬ 
elasticity.  Reference  1  defines  this  new  field  of  hydroelasticity.*  In 
two  of  the  first  stability  studies,*’®  the  classical  translational  and 
rotational  degrees  of  freedom  of  hydrofoils  were  considered.  Henry,  et  al., 
adapted  theory  from  the  field  of  aeroelasticity  and  investigated  the  static 
and  oscillatory  stability  of  an  idealized,  fully  submerged  foil.  They 
concluded  that  flutter  is  very  unlikely  for  the  small  values  of  the  mass- 
density  ratio  1^  typical  of  hydrofoils.  In  a  parallel  study,  however, 
Abramson  and  Chu®  were  more  concerned  about  the  possibility  of  flutter. 


*References  are  listed  on  page  24 
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Both  studies  indicated  that  there  exists  a  minimum  value  of  the  mass-density 
ratio  M,  below  which  flutter  would  not  occur.  The  efforts  of  these  early 
investigators  were  hampered  by  the  lack  of  experimental  flutter  data  with 
which  to  compare  their  theoretical  predictions.  In  the  only  conclusive 
reports  of  flutter  tests'*’^  in  water  available  to  them,  no  instability  had 
been  demonstrated.  They,  therefore,  compared  their  theoretical  results 
with  results  of  experiments®  in  which  light  aeroelastic  models  were 
tested  in  wind  tunnels  using  high-density  gases.  The  comparisons  which 
they  made  to  determine  the  accuracy  of  flutter  theory  revealed  certain 
inconsistencies  which  were  not  completely  reconciled.  Later,  the  results 
of  Hilborne's  study^  became  available.  Hilborne  produced  flutter  of  solid, 
cantilevered,  surface-piercing  foils  in  a  rotating-arm  basin.  He  made 
predictions  on  the  basis  of  an  analysis  which  is  seldom  used  in  the  United 
States.  The  predictions  agreed  with  the  experimental  results,  but  the 
data  reported  were  fragmentary  and  insufficient  to  enable  one  to  make 
comparisons  with  the  classical  theory  of  Theodorsen,  which  is  generally 
used  in  this  country.  It  is  known,  however,  that  for  one  of  Hilborne's 
foils,  the  value  of  the  mass-density  ratio  was  in  the  range  considered  by 
Henry,  et  al.  to  be  typical  of  hydrofoils.  Hilborne  s  foils  were  simple 
flat  plates,  tapered  at  the  leading  and  trailing  edges  and  constructed  of 
solid  steel.  His  results*  have  been  considered  inconclusive  for  three 
reasons : 

1.  Circular  path  of  travel  of  the  foils. 

2.  Possible  interference  from  the  wake  of  a  previous  circuit. 

3.  Attendant  cavitation  at  some  speeds  at  which  flutter  was  obtained. 

In  these  early  stability  studies,  no  consideration  was  given  to  the 
possible  destabilizing  effects  of  a  mass  attached  to  the  oscillating  foil. 
Consideration  of  the  stability  of  such  an  oscillatory  system  was  prompted 
by  the  vibration  problem  encountered  on  USS  FORREST  SHERMAN  (DD  931) .  The 
modes  of  vibration  involved  were  identified  by  full-scale  tests  on  destroyers 


*Gruman  Aircraft  Engineering  Corporation  is  now  conducting  tests  with  a 
linear  towing  path  to  verify  Hilborne's  results. 


2 


t 


of  the  class  as  the  three-noded  horizontal  hull  mode,  coupled  with  rotation 
of  the  rudders. In  this  case,  the  stern  vibrated  in  the  horizontal 
athwartship  direction  and  the  rudders  rotated  as  rigid  bodies  with  respect 
to  the  stern  about  vertical  axes  through  the  rudder  stocks.  The  hull- 
translation  and  rudder-rotation  modes  of  vibration  can  be  considered  to 
correspond,  respectively,  to  the  translational  and  rotational  modes  of  the 
classical  flutter  theory.  An  experimental  study,  conducted  at  the  David 
Taylor  Model  Basin, of  a  mass-foil  system  in  water  did  not  produce 
instability,  but  did  indicate  an  approach  to  flutter.  Jewell^^  discussed 
the  possible  destabilizing  effect  of  an  external  mass  attached  to  an 
oscillating  foil.  He  pointed  out  that  an  additional  mass,  such  as  the 
effective  mass  of  a  ship  hull,  could  increase  the  effective  value  of  the 
mass-density  ratio  enough  to  permit  or  even  promote  flutter.  McGoldrick'"^ 
concluded  that  the  DD  931  had  suffered  a  flutter  condition,  although  the 
problem  was  complicated  by  indications  of  cavitation  on  the  ship's  rudders. 
This  problem  is  very  interesting  and  controversial,  as  evidenced  by 
McGoldrick's  paper the  discussions  of  the  paper,  and  the  author's  reply. 

The  experimental  study  described  herein  is  an  extension  of  the  work 
presented  in  Reference  11.  The  work  was  conducted  at  the  Taylor  Model 
Basin  in  an  effort  to  provide  documented  data  on  flutter  in  water  and  to 
provide  a  firm  basis  for  comparison  with  theoretical  results.  The  experi¬ 
mental  equipment  and  techniques  are  described,  and  the  experimental  flutter 
data  are  presented.  Also  provided  are  comparative  theoretical  data  computed 
on  the  basis  of  a  classical,  two-dimensional,  unsteady  theory.  This 
theory,  originally  derived  by  Theodorsen,^^  was  used  because  it  is  well- 
known  by  experienced  hydroelasticians  and  can  be  easily  acquired  by 
scientists  or  designers  who  are  unacquainted  with  the  fields  This  work 
is  Identified  by  Bureau  of  Ships  Project  Number  S-F013  02  01,  Task  1719. 
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EXPERIMENTAL  EQUIPMENT 


The  flutter  apparatus  was  basically  the  same  as  that  used  previously^^ 
at  the  Taylor  Model  Basin.  The  apparatus  was  mounted  on  a  towing  carriage 
which  j-uns  over  the  high-speed  towing  basin.  The  carriage  is  self-propelled 
and  has  a  top  speed  of  about  50  knots.  The  basin  is  21  feet  wide,  over 
2900  feet  long,  and  is  filled  to  a  minimum  depth  of  10  feet.  Figures  1, 

2,  and  3  show  the  following  pertinent  features  of  the  apparatus: 

1.  The  control  surface  is  a  relatively  rigid,  hollow  hydrofoil 
(la)  with  an  NACA  0015  cross  section,  an  18-inch  chord  length  2b  and  an 
18-inch  span  length  i. 

2.  The  hydrofoil,  vertical  shaft  (lb),  cage  (Ic)  ,  top  disk  (Id), 
and  sliding  counterweights  (2a)  are  all  rigidly  connected  to  form  the 
rotational  element  with  a  vertical  axis  through  the  center  of  the  top  disk, 
the  shaft  axis,  and  the  forward  quarter  chord  of  the  foil. 

3.  The  rotational  degree  of  freedom  is  provided  by  allowing  the 
rotational  assembly  to  turn  with  respect  to  a  strongback  (3a)  to  which 

it  is  connected  by  flexures  (3b).  The  mass  moment  of  inertia  of  the 
rotational  assembly  with  respect  to  the  axis  of  rotation  is  denoted  by  I. 

The  rotational  flexures  provide  an  elastic  restoring  moment  with  spring 

constant  K  .  The  natural  frequency  ou  of  the  uncoupled  rotational  degree 

CL  -  CX 

of  freedom  is  given  by  m  =  Vk  /I. 

a  a 

4.  The  strongback  and  the  entire  rotational  assembly  can 
translate  horizontally  with  respect  to  the  towing  carriage,  to  which  they 
are  connected  by  another  set  of  flexures  (3c)  .  With  this  construction, 
the  mass  m,  which  oscillates  in  horizontal  translation,  includes  the  mass 
of  the  entire  rotational  assembly. 

5.  The  degrees  of  freedom  can  be  locked  independently  or  the 
amplitudes  of  vibration  can  be  limited  by  use  of  adjustable  mechanical 
stops  (le) . 
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6.  The  foil  shaft  passes  through  a  hole  in  an  8-foot  by  16-foot 
rectangular,  horizontal  surface  plate  shown  in  Figure  3  of  Reference  11. 

The  plate  is  used  to  minimize  surface  effects  and  serves  as  an  end  plate 

at  the  upper  end  of  the  hydrofoil.  The'surface  plate  is  submerged  6  inches. 
There  was  approximately  1/16-inch  clearance  between  the  bottom  of  the  sur¬ 
face  plate  and  the  top  of  the  foil. 

7.  A  "take-up"  spring  (If)  can  be  adjusted,  within  limits,  to 
balance  the  steady  lift  force  and  keep  the  shaft  positioned  near  the  center 
of  the  hole  in  the  surface  plate.  The  translational  flexures  and  the 
ajustable  "take-up"  spring  provide  an  elastic  restoring  force  with  spring 

constant  K,  . 

n 

8.  The  location  of  the  center  of  gravity  (c.g.)  of  the  rotational 

element  can  be  adjusted  without  changing  the  mass  m  or  mass  moment  of 

inertia  I  about  the  axis  of  rotation  by  sliding  the  counterweights  around 

the  edge  of  the  top  disk  of  the  rotational  element.  This  provides  a 

convenient  means  of  varying  the  first  moment  of  mass  S  . 

a 

The  basic  apparatus  is  denoted  as  Configuration  A.  Three  other  con¬ 
figurations  —  B,  C,  and  D  —  are  obtained  by  adding  masses:  to  the  top 
disk  for  B,  to  the  strongback  for  C,  and  to  both  for  D.  These  four  inertia 
configurations  are: 

Configuration  A  is  identical  with  that  described  in 
Reference  11  except  that  a  wider  range  of  values  of  the  mass 
unbalance  are  obtainable.  This  was  accomplished  by  removing 
the  forward  damping  coil  (the  black  cylinder  shown  in  Figure  1 
above  If  but  not  in  Figure  2)  and  replacing  it  with  an  equal 
mass  at  the  same  radius  from  the  rotational  axis.  This  new 
mass  is  part  of  the  sliding  counterweights  (item  2a  of  Figure 
2) .  The  small  differences  between  the  values  of  the  spring 
constants  and  given  here  and  in  Reference  11  are  due  to 
more  careful  calibration  in  the  later  studies.  In  this  work, 
m,  S  ,  and  I  are  defined  as  the  structural  inertias,  whereas 
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in  Reference  11  these  quantities  include  the  fluid  added  masses. 
The  altered  values  of  the  uncoupled  frequencies  ojj^  and  (in  air) 
reflect  these  revised  definitions. 


Configuration  B  was  obtained  by  rigidly  attaching  a 

71. 5- lb  weight  (Item  2b  of  Figure  2)  to  the  top  disk  of  the 

apparatus.  This  mass  increased  the  values  of  m  and  I  and 

decreased  the  values  of  tju,  and  o)  .  The  mass  was  positioned  to 

ha 

produce  values  of  S  much  greater  than  those  for  Configuration  A. 

cc 


Configuration  C  was  obtained  by  rigidly  attaching  the  71. 5- lb 

weight  to  the  strongback.  This  increased  the  value  of  m  and 

decreased  the  value  of  from  those  for  Configuration  A  and 

made  them  equal  to  their  values  for  Configuration  B.  Since 

this  mass  was  not  attached  to  the  rotational  element,  it  did 

not  affect  the  values  of  I,  uj  ,  or  S  ,  which  are  exactly  the 

a  a 

same  as  for  Configuration  A. 


Configuration  D  was  obtained  by  attaching  the  71.5-lb 

weight  to  the  top  disk  and  a  51. 5- lb  weight  to  the  strongback 

(Item  2c  of  Figure  2).  This  configuration  is  similar  to 

Configuration  B  in  that  the  values  of  I ,  u)  ,  and  S  are  the 

Cl  Cl 

same.  The  5 1.5- lb  weight  further  increased  the  value  of  m 

and  decreased  the  value  of  u;,  . 

h 

This  oscillatory  system  can  be  considered  as  an  idealized  foil  and 

mass,  connected  by  translational  and  rotational  springs  as  shown  in 

Figure  4.  The  dimensionless  parameter  a  is  the  ratio  of  the  distance 

of  the  axis  of  rotation  aft  of  midchord  to  the  foil  semi-chord  length. 

The  dimensioness  distance  from  the  axis  of  rotation  to  the  center  of 

gravity  (c.g.)  location  is  denoted  by  x  .  It  is  obtained  by  dividing 

Cl 

the  mass  unbalance  by  the  product  of  the  semi-chord  b  and  the  total 

translational  mass  m.  Thus  x  =  S  /mb  and  is  positive  for  c.g.  locations 

a  a 

aft  of  the  axis  of  rotation- 
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Strain  gages  were  bonded  to  the  translational  and  rotational  flexures. 
The  gage  leads  were  wired  to  form  four-arm  bridges,  and  the  amplified  strain 
signals  were  recorded  to  indicate  the  vibratory  displacements.  Carriage 
vibrations  were  measured  with  accelerometers  and  were  recorded.  An  inter¬ 
rupter  device  yielding  one  signal  for  every  10  feet  of  carriage  travel  was 
used  in  conjunction  with  timing  signals  to  determine  the  carriage  speed. 

Limitations  in  the  equipment  complicated  the  experiment.  Two  pertinent 
limitations  were: 

1.  Above  15  knots,  the  transverse  displacement  due  to  the  steady 
lift  force  became  so  large  that  the  foil  shaft  was  pushed  against  the  edge 
of  the  hole  in  the  surface  plate.  Thus,  the  speed  range  was  limited  to 
those  at  which  oscillatory  displacements  of  measurable  amplitude  could  be 
obtained.  Except  for  this  speed  restriction,  the  limited  displacements 
were  not  considered  objectionable.  In  fact,  the  edges  of  the  hole  acted 

as  stops  and  allowed  data  to  be  taken  under  unstable  vibratory  conditions 
without  causing  damage  to  the  equipment. 

2.  At  certain  speeds,  the  frequency  of  vibratory  sources  on  the 
carriage  coincided  with  the  frequency  of  the  apparatus  vibrations.  A 
graph  of  carriage  excite  'on  frequencies  is  shown  in  Figure  5.  Reference 
11  contains  a  detailed  discussion  of  these  vibrations. 

One  other  point  should  be  mentioned  here.  During  a  static  calibration 
it  was  noted  that  the  translational  displacement  at  the  bottom  of  the  foil 
was  slightly  greater  than  the  displacement  at  the  top  of  the  shaft.  This 
difference  is  attributed  to  the  geometry  of  the  translational  flexures 
which  gave  the  effect  of  a  rotational  displacement  about  a  horizontal 
fore  and  aft  axis  located  more  than  20  feet  above  the  foil.  Since  the 
maximum  displacement  at  the  tip  of  the  foil  was  limited  to  about  0.5  in., 
the  angular  motion  could  be  neglected  and  the  foil  motion  considered  to 
be  transverse  translation.  Because  of  the  slightly  rotational  motion  of 
the  foil  in  the  translational  degree  of  freedom,  the  value  of  the  linearized 
spring  constant  IC  (the  ratio  of  a  statically  applied  force  to  the  resulting 
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displacement)  differed  at  the  top  and  bottom  of  the  foil.  The  value  of  K 
was  therefore  determined  by  computing  measured 
frequency  and  m  is  the  mass  obtained  by  weighing. 


The  foil  was  hollow  and  was  filled  with  water  for  both  calibrations 
in  air  and  testing  in  the  basin.  The  increased  damping  due  to  the  water 
in  the  foil  was  negligible,  amounting  to  less  than  0.001  of  critical  damping. 
With  the  above  exceptions,  the  calibration  procedures  were  the  same  as 
those  described  in  Reference  11.  The  values  of  the  dimensional  constants 
and  the  dimensionless  parameters  of  the  apparatus  are  given  in  Tables  1, 

2 ,  and  3 . 


EXPERIMENTAL  RESULTS 


The  specific  objectives  of  the  experiment  were  to: 

1.  Measure  the  oscillatory  damping  over  the  range  of  values 

of  speed  U  and  mass- unbalance  S^  for  as  many  configurations  as  were  feasible. 

2.  Determine  the  quantitative  relation  between  the  reduced 

flutter  speed  b  and  the  dimensionless  distance  from  the  axis  of 

fa 

rotation  to  the  center  of  gravity  x  for  each  value  of  the  mass-density  ratio 

a 

q,  =  m/npb^X  insofar  as  possible. 


A  secondary  objective  was  to  measure  the  frequency  UJ  of  vibration  along 
with  the  damping.  No  attempt  was  made  to  determine  the  ratio  of  the 
amplitudes  of  vibration  of  the  two  degrees  of  freedom  or  to  determine  the 
phase  relation  between  them. 

The  angle  of  attack  of  the  foil  was  preset  at  1.9  degrees.  At  this 
angle,  incipient  cavitation  is  expected  to  occur  at  the  top  of  the  foil 
at  speeds  over  25  knots.  Since  the  speed  was  limited  to  approximately 
15  knots  for  reasons  previously  noted,  it  is  quite  certain  that  the  foil 
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TABLE  1 


Constants  Independent  of  Configuration 


Symbol 

b 

K 

a 

"h 

i 


Value 


Dimensions 


9.00 

3.72  X  10^ 
894 
18.0 


in . 

lb- in . / rad 
lb/ in . 
in . 


TABLE  2 


Parametric  Values  Dependent  on  Configuration 


Symbol 


Value 


Dimens  ions 


A 

B 

C 

D 

m 

1.18 

1.39 

1.39 

1.52 

lb-sec^ / in . 

I 

49.1 

55.3 

49.1 

55.6 

lb-in. -sec® 

(jU 

27.5 

25.9 

27.5 

25.9 

rad/ sec 

a 

27.5 

25.4 

25.4 

24.2 

rad/ sec 

S 

1.49 

2.25 

1.49 

2.25 

lb- sec" 

a 

1.63 

2.40 

1.63 

2.40 

1.77 

2.54 

1.77 

2.54 

1.88 

2.66 

1.88 

2.66 

1.99 

2.77 

1.99 

2.77 

2.16 

2.86 

2.16 

2.93 

2.26 

2.-93 

2.26 

3.04 

3.04 

r 
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TABLE  3 


Dimensionless  Flutter  Parameters 


Symbol 


Value 


a 


-0.5 


A 

B 

C 

D 

ji  =  m/:Tob^i 

2.76 

3.25 

3.25 

3.43 

r  =  ./l/mb^ 
a  ^ 

0.717 

0.703 

0.660 

0.672 

Q  = 

1  a 

1.00 

0.981 

0.924 

0.934 

X  =  S 
a 

i  /mb 
a 

0 . 140 

0.180 

0.119 

0.164 

0.154 

0.192 

0.130 

0.175 

0.167 

0.203 

0.142 

0.186 

0.177 

0.213 

0.150 

0.194 

0.187 

0.221 

0.159 

0.202 

0.203 

0.229 

0.173 

0.214 

0.213 

0.234 

0.181 

0.222 

0.243 
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was  fully  wetted  during  these  tests.  The  adjustable  stop  for  the  transla¬ 
tional  degree  of  freedom  was  set  to  coincide  with  the  maximum  displacement 
permitted  by  the  size  of  the  hole  in  the  surface  plate.  A  short-circuit 
indicator  was  installed  on  this  stop  for  nearly  all  of  the  stability  runs, 
and  the  interrupted  signal,  if  it  occurred,  was  recorded  with  the  other 
data.  There  was  no  indication  that  the  rotational  displacement  was  limited 
by  any  mechanical  stops.  Each  stability  test  was  conducted  with  the 
carriage  moving  at  essentially  constant  speed. 

Sample  records  of  the  vibration  data  are  shown  in  Figure  6.  Figure 
6a  is  a  complete  record  showing  carriage  vibrations,  and  speed  and  time 
signals  as  well  as  unstable  system  oscillations.  The  nature  of  each  of  the 
three  types  of  oscillatory  stability  is  clearly  shown: 

Figure  6b  Stable  oscillation  (positively  damped) . 

Figure  6c  Neutrally  stable  oscillation  (undamped) . 

Figure  6a  Unstable  oscillation  (negatively  damped) . 

The  unstable  vibrations  were  evidently  self-excited  since  they  required 
no  excitation  other  than  the  normal  background  vibrations  of  the  carriage. 
It  was  necessary  to  stop  these  vibrations  externally  and  to  allow  them  to 
increase  in  amplitude  again  in  order  to  measure  the  rate  of  increase  for 
the  determination  of  the  magnitude  of  the  (negative)  damping.  Otherwise, 
the  unstable  oscillations  would  continue  with  the  amplitude  allowed  by 
the  stop  until  the  end  of  the  run  when  the  forward  speed  was  decreased. 

When  the  vibrations  were  stable  or  neutrally  stable  it  was  necessary  to 
excite  them  by  striking  the  apparatus  with  a  timber.  The  amplitude  of 
the  oscillations  would  then  instantaneously  become  greater  than  the 
background  level,  and  then  decrease  exponentially  in  the  case  of  stable 
vibrations  or  remain  constant  in  the  case  of  neutrally  stable  vibrations. 

In  reducing  the  data,  only  those  oscillations  whose  amplitudes  were 
greater  than  the  background  level  and  less  than  those  permitted  by  the 
stop  were  considered.  Any  small  variations  in  the  amplitude  envelope 
were  neglected  by  fairing  through  the  amplitude  envelopes.  The  measured 


11 


logarithmic  damping  decrements  (or  increments)  for  each  degree  of  freedom 
were  converted  to  dimensionless  damping  ratios  in  the  same  manner  as  for 
a  single  degree  of  freedom  system: 


(c/C  w,  =  o— 

c  translation  2rrn 


hoCt) 

h  (t+nr)/ 
o 


or 


(C/C  ) 

c  rotat ion 


2nn 


In 


ao(t)  ^ 
ao(t+nT)  / 


where 


C  is  the  effective  linearized  damping  coefficient 

is  the  effective  critical  damping  coefficient  of  the  foil 
system 

n  is  the  number  of  periods  between  the  measured  amplitudes 

h^(t)  is  the  translational  amplitude  at  time  t 
a^(t)  is  the  rotational  amplitude  at  time  t 
T  is  the  period  of  oscillation 


The  values  of  the  damping  ratios  and  frequencies  of  vibration  o)  are  plotted 

in  Figure  7  as  functions  of  U  for  each  value  of  S  .  All  data  shown  in 

a 

Figure  7  were  obtained  for  oscillations  about  a  preset  angle  of  attack  of 

a  =  1.9  degrees  except  for  Configuration  B  at  S  =  2.25  lb-sec^,  shown 
s  cc 

in  Figure  7c.  In  this  case,  the  preset  angle  of  attack  was  0.6  degree. 

Data  for  which  the  frequency  of  vibration  coincided  with  a  carriage 
excitation  frequency  are  marked  with  flags.  Often  the  effect  of  frequency 
coincidence  was  to  reduce  the  apparent  value  of  the  damping  ratio. 


The  experimental  critical  flutter  speeds  and  frequencies  were 

obtained  as  the  values  for  zero  damping  from  the  damping  curves.  No 

critical  conditions  were  observed  for  Configuration  C  and  only  one  case 

was  observed  for  Configuration  A.  Unfortunately,  for  Configuration  A  at 

S  =2.26  lb-sec^,  the  flutter  frequency  coincided  with  a  carriage 
a 

excitation  frequency.  The  authors  believe,  however,  that  the  indicated 
critical  flutter  speed  is  approximately  correct. 
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For  Configuration  B,  where  S  =  2.25  lb- sec®  ,  the  critical  flutter 

a 

speed  was  obtained  by  extrapolation  of  the  data;  no  actual  instability  was 

observed.  For  Configuration  B,  where  S  =  2.40,  2.54,  and  2.66  lb-sec  , 

cc 

and  for  Configuration  D,  where  S  =2.77  and  2.93  lb-sec  ,  data  were 

obtained  for  only  a  limited  range  of  speeds  near  .  Accurate  interpolation 

of  these  damping  data  to  determine  was  somewhat  difficult.  These 

values  of  are  therefore  considered  to  be  less  accurate  (within  10  percent) 

than  any  other  values.  The  most  accurate  values  of  U_  were  obtained  for 

the  following  values  of  S  :  For  Configuration  B,  S  =  2.77,  2.86,  2.93, 

a  a 

and  3.04  lb-sec®  ,  and  for  Configuration  D,  S  =  3.04  lb-sec®.  The  accuracy 

a 

of  these  values  is  estimated  to  be  within  5  percent. 

The  experimentally  determined  values  of  and  are  listed  in 

Table  4.  The  dimensionless  values  of  the  critical  flutter  speed  will  be 
shown  later  when  the  experimental  data  are  compared  with  the  theoretical 
data . 
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TABLE  4 


Flutter  Speed  and  Frequency  Values 
for  Each  Conf i.qurat  ion 

E^  er  imental 


Critical 


Config. 

lb-  sec^ 
A  2.26 

B  2.25 

2.40 
2.54 
2.66 
2.77 
2.86 
2.93 
3 . 04 

D  2.77 

2.93 
' '  3.04 


in. /sec 

rad / sec 

301 

30.0 

300 

28.2 

252 

28.6 

264 

28.5 

244 

28.8 

253 

28.6 

250 

29.2 

239 

29.4 

229 

29.8 

300 

27.9 

278 

28.1 

274 

28.2 

Theoretical 


in . /sec 

rad/sec 

273 

30.6 

270 

28.0 

256 

28.3 

244 

28.5 

239 

28.7 

235 

28.9 

231 

29.1 

230 

29.2 

227 

29.4 

266 

27.8 

258 

28.0 

254 

28,2 

14 


THEORY 


Tlie  analysis  of  Che  sCability  of  the  oscillatory  system  used  in  the 
experiment  was  based  on  the  classical  two-dimensional  flutter  analysis  of 
aeroelasticity  developed  by  Theodorsen/^  This  analysis  is  appropriate 
for  this  system  since  the  modes  of  vibration  are  very  well  defined  and  the 
foil  oscillates  as  a  rigid  body  with  elastic  restraints.  The  use  of  the 
two-dimensional  expressions  for  foils  of  small  aspect  ratio  is  not  new  in 
this  work.  For  example,  Widmayer,  et  al.^^  corapardd  their  experimental 
flutter  data  with  reference  theoretical  values  obtained  from  two-dimen¬ 
sional  theory.  Their  models  had  aspect  ratios  ranging  from  2  to  13.  For 
aspect  ratios  from  2  to  6,  the  two-dimensional  theory  predicted  lower  flut¬ 
ter  speeds  than  those  obtained  experimentally.  The  difference  increased 
as  the  aspect  ratio  decreased.  Admittedly,  the  two-dimensional  flutter 
theory  can  be  somewhat  inaccurate  when  used  for  lov;  aspect  ratio  foils; 
however,  it  should  be  recognized  that  a  slight  underestimation  of  the 
flutter  speed  is  desirable.  Abramson  used  two-dimensional  hydrodynamic 
expressions  in  a  study  of  the  stability  of  a  control  surface  on  a  subma- 
rine.  Woolston  and  Castile  also  used  such  expressions  in  their  study. 
Although  they  dealt  with  aspect  ratios  from  4  to  8,  they  showed  that  the 
use  of  the  three-dimensional  expressions  produced  only  slightly  different 
predictions  of  the  critical  speeds.  Reference  17  stated  that  three-dimen¬ 
sional  effects  are  less  significant  at  high  values  of  the  reduced  frequen¬ 
cy,  e.g.,  those  obtained  in  the  present  case, than  at  low  values.  Due  to 
the  complexity  of  the  three-dimensional  expressions,  their  use  often  re¬ 
sults  in  an  uneconomical  expenditure  of  effort. 

The  equations  of  motion  for  the  system  shown  in  Figure  4  can  be 
written:* 


mh  -f  S^a  +  Kj^h  =  -  L  [1] 

S  h  +  la  +  K  a  =  M 
a  a 

*Dots  indicate  differentiation  with  respect  to  time. 
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where 


h  is  the  translational  displacement  of  the  axis  of  rotation 
from  the  foil  equilibrium  position, 
a  is  the  angular  displacement  of  the  foil  from  its  equilib¬ 
rium  position,  and 
L  and  M  are  the  unsteady  hydrodynamic  lift  force  and  moment  refer¬ 
red  to  the  axis  of  rotation. 


Structural  damping  forces  have  been  neglected  in  the  derivation  of  these 
equations.  A  brief  survey  of  aeroelastic  flutter  calculations  indicates 
that  inclusion  of  such  forces  is  seldom  necessary  .  For  the  subject  ap¬ 
paratus,  at  rest,  the  damping  ratio  C/c^  in  water  is  less  than  0.01  for 
each  uncoupled  degree  of  freedom. 


If  simple  harmonic  motion  is  assumed, 
Equations  [1]  and  [2]  can  be  written: 


a  =  a  e 
o 


ioJt 


5 


and 


-moP  h  -  S  tip  a  +  mtiP  h  =  -  L 
a  h 


[3] 


-S  tiPh  -  luPa  4-  luPa  =  M  [4] 

a  a 

The  expressions  used  to  represent  the  oscillatory  lift  force  L  and  moment 
M  are  those  applicable  to  two-dimensional  foils.  The  foil  used  in  the  ex¬ 
periment  was  of  unit  aspect  ratio;  however,  with  the  surface  plate  at  one 
end,  the  effective  aspect  ratio  was  very  nearly  two.  It  is  intended  to 
demonstrate  the  validity  of  the  use  of  the  two-dimensional  expressions  a 
posteriori  when  the  theoretical  data  are  compared  with  the  experimental 
data , 


The  two-dimensional,  unsteady  lift  and  moment  expressions  given  by 
Theodorsen^'^were  multiplied  by  the  length  of  the  foil  I  to  obtain  the 
total  unsteady  hydrodynamic  loading  expressions 


L  =  npb^  ^[ii  +  Ua  -  £ba]  +  2TTpbiUC(k)  [h  +  Ua  +  (^  -  a)ba] 
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[6] 


M  =  npb^i[ah  -  -  a)Ua  -  (1/8  +a®)ba] 

+  2npb^^UC(k)(a  +  ^)[h  +  Ua  +  (^  -  a)ba] 


which  can  also  be  written  as 


L  =  -  npb^^oi?  fL.  h/b  +  [L  -  L.  (^  +  a]a1 
n  a  n  — 


[7] 


M  =  npb^^u?{[M^-  +  a)]h/b 

+  [M^  -  (L^  +  +  a)  +  +  a)2]a} 


[8] 


where  h  and  a  have  been  assumed  to  vary  sinusoidally  with  respect  to  time. 

In  these  last  expressions,  Equations  [7]  and  [8],  the  notation  of  Smilg 
1  8 

and  Wasserman  was  used: 


=  1  -  2iC(k)/k 

^  -  [1  +  2C(K)]i/k  -  2C(k)/k^ 

Mj,  =  i  [91 

=  3/8  -  i/k 
cx 


Here 


k  = 


oib 


U 


and  C(k)  = 


(s) 

Hr^k) 


Hi‘^’(k)  +  inf  (k) 


[10] 


where  and  H  ^  are  Hankel  functions  of  the  second  kind,  the  first 

1  o 
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and  zero  orders,  respectively. 


Since  the  axis  of  rotation  of  the  TMB  flutter  apparatus  is  located 

at  the  forward  quarter-chord,  a  =  -  ^  and  several  of  the  hydrodynamic 

terms  in  Equations  [7]  and  [8]  are  zero.  Substitution  of  these  simplified 

expressions  into  Equations  [3]  and  [4]  and  division  of  Equation  [3]  by 

TTpb^i'O^e^  and  Equation  [4]  by  e^^*",  gives  the  following  flutter 

equations  in  terrnis  of  the  dimensionless  parameters  p,  0,  and  r^  and 

the  dimensionless  amplitudes  h  /b  and  a  : 

0  o 

f  q  [  1  -  C?  (u)  /u))^  ]  +  L,  Ih  /b+^xq+L^Ct.  =  0 
^  ^  ^  a  ^  '  ho  a  a  0 

^)ix  +  11  ^h  /b  +  ^.ir  ^[1  -  (aj  /ou)2  ]  +  M  ]a  =0 
a  ho  a  a  a  o 


[11] 

[12] 


These  equations  express  relations  between  the  various  parameters  at 

the  critical  flutter  speed  and  frequency.  To  investigate  the  stability  of 

the  system  under  conditions  other  than  critical,  a  damping  coefficient  g 

was  used.  The  quantity  (o)  /o))^  in  Equations  [11]  and  [12]  is  replaced  by 

GC 

the  dimensionless  quantity 


Z  =  Cua/wf  (1  +  ig)  [13] 

where  g  can  be  interpreted  as  representing  the  amount  of  damping  required 
to  produce  neutral  stability. 


The  homogeneous  simultaneous  equations,  Equations  [11]  and  [12],  have 

nontrivial  solutions  only  when  the  determinant  of  the  coefficients  of  h  /b 

0 

and  a  is  zero: 
o 

-  [|i.Cf  (nr^  +  M  )  +  ur^  (M,  +  L,  )]Z 


-  [|.LX^  +  K^]  [M-x^  +  L^]  +  [p,  +  Lj^]  [M-r^^  +  M^]  =  0 


[14] 
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This  quadratic  can  be  solved  for  the  complex  quantity  Z  and,  hence,  the 
damping  coefficient  g  and  the  frequency  UJ.  in  the  solution  of  this  com¬ 
plex  quadratic,  values  of  •■H  and  g  are  computed  for  assumed  values  of  k. 

The  corresponding  speeds  are  then  found  from  U  =  UJb/k. 

The  values  of  g  and  tu  are  plotted  as  functions  of  U  for  each  set  of 
values  of  [j,,r^,x^,  and  Q  in  Figure  8.  The  value  of  g  is  negative  for 
stable  oscillations  and  positive  for  unstable  oscillations.  As  usual  in 
flutter  calculations,  only  one  of  two  roots  of  the  complex  quadratic  in¬ 
dicates  an  instability.  Of  the  two  roots,  the  values  of  g  and  <X)  shown  in 
Figure  8  correspond  to  the  mode  of  vibration  which  is  unstable  for  suffi¬ 
ciently  large  values  of  x  .  The  theoretical  values  of  U..  and  were  ob- 

a  f  f 

tained  by  interpolation  (corresponding  to  g  =  0)  in  the  same  manner  as  the 
experimental  values. 


DISCUSSION 

Although  the  damping  parameters  g  and  C/C^  are  related,  they  are  dif¬ 
ferent  in  concept;  therefore  no  quantitative  comparison  of  their  values  is 
attempted.  It  is  noted,  however,  that  the  corresponding  damping  versus 
speed  curves  are  qualitatively  similar.  The  reader  can  compare  values  of 
the  experimental  and  theoretical  frequencies,  but  this  comparison  is  con¬ 
sidered  to  be  of  secondary  interest  at  present.  Of  utmost  importance  here 
is  a  quantitative  comparison  of  the  experimental  and  theoretical  values  of 
the  flutter  speed.  For  this  purpose,  is  shown  plotted  against  S  for 
each  configuration  in  Figure  9.  Table  4  also  lists  and  for  each  con¬ 
figuration.  In  addition,  the  reduced  flutter  speed  parameter  ^  is 

shown  as  a  function  of  the  dimensionless  distance  from  the  axis  of  rotation 
to  the  center  of  gravity  x  in  Figure  10.  It  is  noteworthy  that  flutter 
occurred  at  the  relatively  low  speeds  of  11.5  to  16  knots. 

It  is  clear  from  Figures  9  and  10  that  an  increased  value  of  or 
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X  results  is  a  decreased  critical  flutter  speed  for  the  range  of  values 

Ou 

of  system  parameters  (u,,r  ,0)  considered.  This  destabilizing  effect  of 

cx 

•  increasing  S  is  common  knowledge  among  aeroelasticians .  The  analytical 

results  agree  with  the  experimental  results  in  showing  this  trend.  The 
critical  flutter  speed  can  be  larger  for  much  larger  values  of  and  can 
remain  finite  for  negative  values  of  S^.  The  predicted  critical  flutter 
.  speeds  are  either  conservative  (a  maximum  of  12  percent ‘less  than  the  ex¬ 
perimental  critical  speed)  or  lie  within  the  experimental  error.  In  fact, 
only  three  of  the  theoretical  critical  flutter  speeds  differ  from  the  ex¬ 
perimental  speeds  by  more  than  the  estimated  experimental  accuracy.  The 
authors  believe  that  the  results  of  the  comparison  just  stated  definitely 
justify  the  use  of  the  unsteady  two-dimensional  flutter  theory  for  the 
values  of  p,,  a,  and  k  involved  in  this  study  to  determine  critical  flutter 
speeds  of  similar  configurations. 

In  these  results,  the  reduced  frequency  k  at  the  critical  flutter 

speed  and  frequency  (t^J^b/U^)  ranges  from  a  theoretical  value  of  0.833  for 

Configuration  B,  where  S  =  2.54  lb-sec^,  to  an  experimental  value  of  1,17 

cx 

for  Configuration  B,  where  S  =  3.04  lb-sec^  .  These  are  very  high  values 

CX 

with  respect  to  those  prevalent  in  aeroelasticity  where  typical  values  of 
k  are  of  the  order  of  0.1.  It  was  reported  that  several  two-dimensional 
aerodynamic  terms  in  Theodorsen's  expressions  for  the  unsteady  lift  and 
moment  differ  from  experimental  values  for  k  above  0.4.  Therefore,  one 
might  expect  poor  agreement  between  theoretical  and  experimental  flutter 
results.  In  the  current  study,  however,  the  theoretical  predictions 
shov/n  in  Figures  9  and  10  are  in  good  agreement  with  the  experimental 
flutter  speeds  and  frequencies.  These  results  appear  to  reinforce  the 
statement  made  previously  that  three  dimensional  effects  are  less  signif¬ 
icant  at  high  values  of  the  reduced  frequency  than  at  low  values. 

Since  both  r  and  0  were  changed  with  p  in  this  experiment,  no  con- 
cc 

elusions  concerning  the  effect  of  p  on  the  trends  could  be  obtained. 
Although  the  predicted  and  experimental  speeds  agree,  there  remains  the 
possibility  that  the  flutter  theory  used  herein  would  yiel’d  unconserva¬ 
tive  critical  speed  predictions  for  lower  values  of  p.  This  would 
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correspond  to  the  inaccuracy  of  flutter  predictions  of  Woolston  and 

Castile  for  values  of  u,  near  Values  of  u,  were  computed  for 

cr 

two  values  of  x  for  Configuration  B.  For  x  =  0.192,  u,  =  2.44  and  for 
a  a  cr 

=  0.243,  =  1.93,  as  shown  in  Figure  11.  The  experimental  value  for 

p  (3.25)  was  too  much  larger  than  to  indicate  the  accuracy  of  the 

theoretical  predicticnof  near  u  .  Due  to  the  limitations  of  the  exist- 

f  cr 

ing  flutter  apparatus,  the  authors  could  not  investigate  this  matter.  Mr. 
C.  J.  Henry  of  the  Davidson  Laboratory  is  investigating  this  dependence 
on  p  both  by  experiment  in  water  and  by  analysis. 

It  appears  appropriate  to  consider  what  this  and  earlier  studies  in¬ 
dicate  concerning  the  possibility  of  flutter  occurrence  on  displacement 
vessels  and  hydrofoil  craft.  The  control  surfaces  of  ships,  such  as  the 
DD  931  rudders,  are  usually  mass  balanced;  that  is,  x  is  nearly  equal  to 
zero.  Previous  studies  show  that  flutter  occurs  at  relatively  high  values 
of  the  reduced  flutter  speed  U^/u^^^'^for  x^  =  0;  see  Figure  12,  taken  from 
Reference  19.  In  addition,  the  flutter  speed  is  relatively  high  for  large 
values  of  p,  such  as  p  =  36,  which  is  the  effective  value  calculated  for 
DD  931.  This  indicates  that  until  speeds  of  surface  ships  or  submarines 
are  increased  significantly,  they  are  not  in  danger  of  hull-control-sur¬ 
face  flutter,  especially  where  no  free-surface  or  cavitation  effects  are 
involved.  In  spite  of  this,  both  the  experimental  results  and  Figures 
12  d,  e,  and  f  show  that  flutter  could  occur  at  low  speed  on  a  configura¬ 
tion  where  the  center  of  gravity  of  a  control  surface  lies  somewhat  aft 
of  the  axis  of  rotation.  In  this  case,  the  minimum  flutter  speed  occurs 
for  a  value  of  near  unity,  which  is  the  value  derived  for  the  DD  931. 

When  the  possibility  of  hydrofoil-craft  flutter  is  considered,  the 
indication  is  more  serious.  Although  the  value  of  p  is  slightly  below 
p  ^  for  ordinary  foil  sections,  the  inclusion  of  propeller  nacelles  or 
large  end  plates  on  hydrofoil  configurations  could  effectively  increase 
p  to  values  above  p^^  In  this  case,  hydrofoil  craft  could  be  susceptible 
to  flutter  because  of  their  high  speed  and  because  of  the  many  possible 
coupled  modes  of  vibration  of  the  flexible  foils,  struts,  control  surfaces, 
and  hull.  Using  aeroelastic  experience  with  troublesome  "Tee  Tail" 
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flutter  problems  as  a  guide,  the  combination  of  a  horizontal  foil  attached 
to  a  vertical  strut  would  be  particularly  susceptible  to  flutter.  In  any 
case,  if  flutter  occurs  at  a  value  of  p,  slightly  greater  than  it  will 

occur  near  the  least  possible  critical  speed,  as  shown  in  Figure  12. 

The  pertinent  points  can  be  summarized  as  follows; 

1.  Flutter  can  occur  on  displacement  vessels  at  high  speeds  in  the 
coupled  hull-bending,  control-surface-rotation  mode  of  vibration  for  suf¬ 
ficiently  large  values  of  x^.  The  critical  flutter  speed  is  reduced  for 
values  of  Q  near  unity 

2.  On  hydrofoil  craft,  if  coupling  occurs  between  any  of  the  many 
possible  vibratory  modes  and  if  the  effective  value  of  p  is  greater  than 
p,^^,  then  flutter  is  quite  probable. 

It  should  be  emphasized  that  these  experiments  utilize  idealized 
systems  representing  specified  modes  and  couplings.  Other  vibratory  modes 
and  couplings  must  be  studied,  together  with  the  special  effects  of  the 
free  surface,  cavitation,  and  waves.  These  last  three  phenomena  introduce 
new  boundary  conditions  which  distinguish  hydroelastic  from  aeroelastic 
phenomena.  Analytical  studies  of  these  effects  indicate  that  cavita¬ 
tion  will  seriously  worsen  the  flutter  problem.  The  only  experimental 

g 

work  on  these  effects  known  now  is  that  of  Hilborne.  Until  the  influences 
of  these  special  effects  are  determined  experimentally,  the  seriousness 
of  the  threat  of  flutter  to  naval  configurations  will  remain  in  doubt.  It 
is  therefore  recommended  that  experimental  studies  be  pursued  to  cover  the 
wide  range  of  current  and  projected  hydrofoil  operating  conditions. 
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CONCLUSIONS 


Flutter  of  a  control  surface  has  been  definitely  demonstrated  in 
water.  Although  flutter  is  generally  considered  to  be  a  high-speed 
phenomenon,  it  was  produced  at  the  relatively  low  speed  of  11.5  knots. 

A  simple  flutter  theory  was  shown  to  be  adequate  to  predict  the  phenomenon 
in  this  study  but  some  doubt  remains  about  its  accuracy  for  other  struc¬ 
tural  configurations  and  flow  conditions. 

These  results  indicate  that  flutter  could  be  an  extremely  serious 
hydro-elastic  problem,  particularly  for  hydrofoil  craft. 
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Figure  1  -  TMB  Flutter  Apparatus 
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Figure  3  -  Sectional  Elevation  of  TMB  Flutter  Aonaratus 


Figure  fi  -  Sample  Oscillatory  Data 


32 


Figure  6b  —  Stable  Oscillation 


Figure  6c  —  Neutrally  Stable  Oscillation 
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Fi^uro  10  -  nimensionlcss  Critical  Flutter  Parameters 
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Figure  12  —  Theoretical  Critical  Flutter  Speeds 


INITIAL  DISTRIBUTION 


Copies 

8  CHBUSHIPS 

3  Tech  Info  Br  (Code  335) 

1  Prelim  Des  (Code  420) 

1  Lab  Mgl (Code  320) 

3  Hull  Des  (Code  440) 

1  CHBUWEPS 

(Code  RAAD-222) 

Alin:  Mr.  D.  Michel 

2  CHONR 

I  Flu  Dyn  Br  (Code  438) 

I  CDR,  USNOL,  While  Oak,  Md 

1  CDR,  USNOTS,  China  Lake 

2  COR,  USNOTS,  Pasadena  Anne*  (Code  P-508) 

1  CDR,  ARESOEVCOM 

Alin:  Mech  Br,  AF  Ollicc  of 
Sti  Res 

1  CDR,  W-PADEVDIV,  Aircrall  Lah 
Alin:  Mr.  W.  Mykylow,  Dyn  Br 

10  COR,  ASTIA 

1  CHAFSWP 

1  SNAME,  N.Y. 

1  NRC,  Canada 

1  NRC,  Washinglon 

1  ADMIN,  U.S.  Maiilime  Adm 

Alin:  Ml.  R.P,  Godwin,  Acling  CHF, 

OK  of  Res  S  Dev 

2  DIR,  Nall  BuSland 

1  Dr.  G.B.  Schubauei,  CHF,  Flu 
Mech  Sec 

1  Or.  J.M,  Fiankland,  ConsullanI 

2  DIR,  Langley  Res  Cli,  NASA,  Hydi  Div 

1  Ml.  I.E.  Garrick 
1  Ml,  O.J.  Marlon 

3  DIR,  DL,  SIT,  Hoboken 

1  Ml.  C.J,  Henry 
1  Ml.  S.  Tsakonas 

2  DIR,  ORL  Penn  Slalc 

1  Dr.  M.  Sevik 

3  DIR,  Depl  of  Mech  Sci,  SIV  Res  InsI 

San  Anlonio 
1  Dr.  H.N.  Abramson 
1  Ml.  G.  Ransleben 

2  DIR,  Iowa  InsI  ol  Hydiau  Res,  Slalc  Univ 

of  Iowa 

1  Prof.  L.  Laiidwebei 

3  DIR,  SI  Anlhony  Falls  Hydiau  Lab,  Univ 

ol  Minn,  Minneapolis 
1  Prof.  B.  Silbeiman 
1  Ml.  J.N.  WcUfI 

1  DIR,  AppI  Phys  Lab,  Johns  Hopkins  Univ 
Silver  Spring 

1  Elec  Boal  Div,  Genl  Dyn  Coip,  Giolon 

Alin:  Ml.  Robeil  McCandliss 

2  Slanlord  Univ,  Depl  ol  Malh,  Slanfoid 

I  Or.  B.  Perry 
1  Dr.  E.Y.  Hsu 


Copies 

1  Univ  of  Calif,  Depl  ol  Engin,  InsI  of  Engin 
Res,  Berkeley,  Alin:  Dr.  J.V.  Wchausen 

3  Calif  InsI  ol  Tech,  Pasadena 
1  Or.  M.S.  Plessel 
I  Dr.  T.Y.  Wu 
1  Dr.  A.J.  Acosla 

3  Mass  InsI  ol  Tech,  Flu  Dyn  Res  Lab, 

Cambridge 
I  Prof.  H.  Ashley 
I  Prof.  M.  Landahl 
I  Prof.  J.  Dugundji 

1  MidwesI  Res  InsI,  Kansas  City 
Alin:  Ml.  Zcydel 

4  Univ  ol  Mich,  Ann  Arbor,  Depl  of  Engin  Mech 

1  Prof.  Jesse  Oimondioyd 
1  P.ol.  C.S.  Yih 
Depl  of  NAME 

I  Piol.  R.B.  Couch,  Dll  E*p  Nav  Tank 
I  Prof.  F.  Michaelson 

1  Piol.  H.A.  Schadc,  Hd,  Depl  ol  Nav  Arch 
Univ  ol  Calil,  Bcikeley 

1  Prof.  A.G.  Sliandhagen,  Depl  ol  Engin  Mech 
Univ  ol  Noire  Dame 

1  Prof.  W.R.  Sears,  Grad  School  ol  Aero  Engin, 
Cornell  Univ,  llhaca 

I  Prof.  J.S.  McNown,  Dean  ol  Engin,  Univ  ol 
Kansas,  Lawrence 

I  Piol.  F.A,  Bibcislein,  Hd,  Depl  ol  Civil  Engin 
Calholic  Univ,  Washinglon 

1  Piol.  B.  Morrill,  Chairman.  Mech  Engin  Depl 
Swailhinoie  College,  Swarlhmoie 

1  Boeing  Airplane  Co,  Seallle 

Alin:  Mr,  M,J.  Turner 

2  Convaii,  San  Diego 

1  Ml.  A.D.  MacLellan,  Sys  Dyn  Gp 
1  Ml.  H.T.  Brooke,  Hydro  Gp 

2  Tech  Res  Gp,  Inc,  Syossel,  Long  Island 
1  Dr.  J,  Kolik 
1  Piol.  S.  Karp 

1  The  Rand  Corp,  Sania  Monica 

Alin:  Or.  B.  Parkin 

2  Cornell  Aero  Lab,  Bullalo 

1  Dr.  Slallei 
1  Ml.  R,  While 

I  Genl  AppI  Sci  Lab,  Inc,  Weslbuiy,  Long  Island 
Alin:  Dr.  F.  Lane 

1  Gibbs  and  Co*.  Inc,  N.Y, 

2  Grumman  Aircrall  Engin  Corp,  Belhpage, 

Long  Island 
I  Ml.  E.  Band 
1  Ml.  C.  Siiulrcs 

1  Grumman  Aircrall  Engin  Corp.  Dyn  Dev  Div 

Babylon 

2  Pres,  Hydronaulics,  Inc,  Rockville 

1  Lockheed  Aircrall  Corp,  Missiles  8  Space 
Div,  Palo  Alio 
Alin:  Mr.  R.W.  Kcrniecn 


Copies 

1  Dr.  P.  Kaplan,  Pres,  Oceanics,  Inc. 

New  York 

4  Natl  Phys  Lab,  Felthanr,  Middlesex,  England 
1  Dir,  Ship  Div 
1  Hd,  Aerodyn  Div 
1  Ml.  A.  Silverleaf 
1  Ml,  D.V.  Hilborne 

2  DIR,  Royal  Aircrall  Estab 

Farnborough,  Hants,  England 
Alin:  Mr.  M.O.W.  Wolle 

1  Dr.  R.  Timman,  Natl  Luchivoortlaboratorium 
Slatcrwag  145,  Amsleidam,  The  Netherlands 
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